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Abstract. Most previous works study the evolution of cooperation in a structured 
^-j" population by commonly employing an isolated single network. However, realistic 

networks are composed of many interdependent ones coupled with each other, rather 
than the isolated single network. In this paper, we consider a network including two 
interdependent networks A and B with the same size, and each node in network A 
probabilistically connects the corresponding node in network B. We introduce the 
public goods game into such network, and study how the probabilistic interconnection 
influences the evolution of cooperation. Simulation results show that there exists an 
intermediate region of interconnection probability leading to the maximum cooperation 
level in the whole network. Interestingly, we find that at the optimal interconnection 
probability the fraction of internal links between cooperators in networks A and B 
is maximal. Also, even if initially there are no cooperators in one interdependent 
network, cooperation can still be promoted by probabilistic interconnection, and the 
cooperation levels in networks A and B can more easily reach an agreement at the 
intermediate interconnection probability. Our results may be helpful in understanding 
the cooperative behavior in some realistic interdependent networks and thus highlight 
the importance of probabilistic interconnection on the evolution of cooperation. 



PACS numbers: 87.23. Ge, 87.23.Kg, 89.75.Fb 
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1. Introduction 

The emergence of cooperation among selfish individuals in real world is still a challenging 
problem in social and biological systems pp. Evolutionary game theory has provided 
a uniform framework to improve our understanding of the emergence and sustenance 
of cooperation among unrelated individuals [21 El H]. The public goods game (PGG), 
as one of the most famous paradigms, is often used for discussing the conflict between 
individuals and common interests [5]. 

In a typical PGG, cooperators (C) contribute an amount c to the public good 
and defectors (D) do not contribute. The total contribution is multiplied by an 
enhancement factor r (r > 1) and then equally distributed by all the members in 
the group. Therefore, Ds obtain the same benefit of Cs with no cost, which confronts 
the individuals with the temptation to defect by taking advantage of the public good 
without contribution. Correspondingly, the Tragedy of the Commons is induced. In 
order to elucidate why cooperators thrive under the exploitation of defectors in PGG, 
many mechanisms have been proposed, such as voluntary participation [61 [7J, social 
diversity El EES E] , punishment pj ECU EGA EES] , migration [EH EE?], reward PS EG], 
coordinated investments [20], the Matthew effect [21], adaptive and bounded investment 
returns [22J, and conditional strategies [23]. In particular, the evolution of cooperation 
in the PGG on complex networks receives more attention recently. It is shown that the 
mesoscale structure plays an important role in promoting the evolution of cooperation 
on complex networks [2J1 125] . 

It is worth pointing out that all the previous works enumerated above are based 
on the limited case of an isolated single complex network, e.g., square lattice, regular 
ring, or scale-free network. However, empirical evidences show that real world networks 
are not isolated or disintegrated, but constructed by lots of interdependent networks, 
which connect and influence one another directly or indirectly [26]. In other words, the 
real world in which we are living is a huge network of networks [271 I2H1 123 EQl [3Tj . 
Especially in social systems, there exists a diverse kinds of related social networks, such 
as the collaboration network and the social network. It is known that, because of the 
social attribute of human beings, these two networks are directly linked together to 
some degree for scientists [32] . Each individual on one network may probably link with 
the corresponding individual on the other network directly. In this situation, the two 
different networks may entangle with each other to different degrees. In particular, 
the relationship between these two interdependent networks can be described by the 
probabilistic interconnection between them. Thus, it is interesting to study the social 
dynamics on the realistic network constructed by at least two interdependent networks 
with probabilistic interconnection, and far less attention has been paid on the evolution 
of cooperation in this type of networks. Moreover, we would like to distinguish the 
interdependent networks from community structure network and hierarchical network 
by network structure and social implication. From the network structure, in a 
community structure network, the connections are tight in communities but sparse 
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Figure 1. Interdependent networks composed of networks A (denoted by red nodes) 
and B (denoted by blue nodes) . Individuals are arranged on two square lattices of the 
same population size. Under the probabilistic interconnection, a fraction of nodes in 
network A are connected with the corresponding nodes in network B. Here, a node 
from one network can connect at most one node from the other network. 



between them [33| 134"] . In a hierarchical network, small groups of nodes organize in 
a hierarchical manner into increasingly large groups, while maintaining a scale-free 
topology [35]. However, in interdependent networks, the interacting layers are connected 
through intra-layer links, which are independent of the interacting networks and can be 
tight or sparse in different situations. The probabilistic interconnection between the 
interconnected networks determines the coupling effect between them |36j. From the 
social implication, the communities in a community structure network probably share 
common properties and play similar roles [37]. They can interact with the others to 
complete one overall functionality. While in a hierarchical network, the individual in 
the center of bigger cluster means being at higher level in the hierarchy [38]. However, 
the networks composing interdependent networks may have different properties or roles. 
These networks interact with each other, influence each other and complement each 
other [39]. Remarkably, recently Wang et al studied the impact of biased utility 
functions on interdependent networks which were connected by utility functions [40J. 
They showed that the benefits of enhanced public cooperation on the two interdependent 
networks are as biased as the utility functions. They emphasized that the positive effect 
of biased utility functions is due to the suppressed feedbacks of individual success, 
which leads to a spontaneous separation of characteristic time scales of the evolutionary 
process on the two interdependent networks. However, we would like to point out that 
interdependent networks in their work are only coupled by the utility functions, and 
indeed they may be also connected by probabilistic interconnection. 

In this paper, we develop a simple model of two interacting networks connected 
through probabilistic interconnection to study the evolution of cooperation on them 
and the coupling effect between them. For simplicity, we employ two identical spatial 
structures A and B of the same size. Moreover, we assume that each node in 
network A connects the corresponding node in network B with probability p. Here 
the interconnection probability p, which controls the number of links between network 
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A and network B, represents the integration of two interconnected networks. For p = 0, 
there are no internal links between the two networks. In other words, networks A and 
B are totally separated. In the opposite limit, that is, p = 1, all the nodes in networks 
A and B are completely connected in order. For < p < 1, the actual number of 
internal links between A and B is subject to a binomial distribution. We find that there 
exists an optimal intermediate region of the interconnection probability p maximizing 
the cooperation level in the whole network. Surprisingly, even if initially there are no 
cooperators in one interdependent network, by means of the coupling effect between 
them, the cooperation level on both networks can more easily reach an agreement at an 
intermediate interconnection probability. 

2. Model 

We investigate the evolutionary PGG on interdependent networks A and B, as shown in 
figure 1. In order to focus explicitly on the impact of the probabilistic interconnection 
and easily compare our results with previous works, the two interdependent networks 
we employed are both L x L square lattices with periodic boundary conditions and 
von Neumann neighborhoods. They have the same population size, Na = Nb- Each 
node in network A connects the corresponding node in network B with probability p. 
Meanwhile, we assume that a node from one network connects to no more than one 
node from the other network. 

Initially, individuals in each network are designated either as a cooperator or a 
defector with equal probability. Depending on the local links between individuals and 
the internal links between interdependent networks, individuals on both networks not 
only engage in five local PGGs which are centered on himself and the nearest neighbors 
on the same network, may but also engage in one long-range PGG which is centered 
on the corresponding node on the other network. In addition, cooperators contribute 
c = 1 to every PGG involved, and defectors contribute nothing. The total contribution 
is subsequently multiplied by the enhancement factor r, and then shared equally by all 
of the group members irrespective of their strategies. The results obtained below are 
from the renormalized PGG enhancement factor n = r/(M + l+p) j6j [9], [201 E]> where 
M = 4 is the number of local nearest neighbors of focal individuals. The payoff of 
defectors engaging in one PGG is Pd = r ■ n c j [M + l+p) = n-n c , and the corresponding 
payoff of cooperators is then P c = P^ — 1, where n c is the number of cooperators in the 
group. After engaging in all the groups, player x is allowed to learn from a randomly 
selected neighbor y including the long-range one. To be specific, player x adopts the 
strategy of the random neighbor y with a probability determined by the difference of 
their payoffs 

W ^ y) = l + eM(P x -Py)/^Y (1) 
where k denotes the noise effect in the strategy adoption process [12] . Following previous 
study [13], we simply set k = 0.5 in this work, and mainly focus on the impact of 
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Figure 2. The frequencies of cooperators among interdependent networks as a 
function of normalized enhancement factor r\ with different values of p in (a) and 
as a function of interconnection probability p with different values of rj in (b) . 



probabilistic interconnection. 



3. Simulation and analysis 

In the following, we show the simulation results carried out on the interdependent 
networks both of size 100 x 100. It is known that network size can strongly influence 
the dynamics of such system [Sj. Thus, we have also used different larger system size 
in our study, e.g., 200 x 200. The simulation indicates that our results are robust 
in larger systems. Initially, cooperators and defectors are randomly distributed on 
both interdependent networks with equal probability. We define that p ac denotes 
the frequency of cooperators of network A, and pb c the frequency of cooperators of 
network B correspondingly. Meanwhile, we define that p c denotes the frequency of 
cooperators in the whole network. In our work, we adopt the synchronous Monte Carlo 
simulation procedure to update the strategies of players. Unless otherwise stated, all 
the simulation results shown below are required up to 3 x 10 4 generations and then 
sampled by another 10 3 generations. The results are averaged over 30 realizations of 
different initial conditions. 

We first plot the frequency of cooperators in the whole system as a function of the 
renormalized enhancement factor n for different values of p in figure 2(a). For p = 0.5, 
we find that the cooperation level is enhanced obviously on a large scale of n and the full 
cooperation state is achieved at n = 0.92. While p = 0, we note that the full cooperation 
state is achieved at n = 1.04, which falls behind p = 0.5. For p = 0.9, at the same values 
of r), the cooperation level is higher than that of p = but lower than that of p = 0.5, 
and the full cooperation state is achieved at n = 0.98. We note that, as shown in figure 
2(a), the probabilistic interconnection between two interacting networks can significantly 
influence the evolution of cooperation in the whole system. In order to study the role 
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Figure 3. The complete p—r/ phase diagram for PGG on the interdependent networks. 
The upper (lower) boundary is the extinction thresholds of cooperators (defectors) 
correspondingly. 

of probabilistic interconnection in promoting cooperation intensively, we investigate the 
frequencies of cooperators as a function of p with different i] in figure 2(b). We observe 
that, for each fixed rj, there exists an optimal intermediate region of p maximizing 
the cooperation level in the whole system. The probabilistic interconnection promoting 
cooperation resembles an interesting resonancelike phenomenon reflected by the optimal 
cooperation level at the intermediate interconnection probability. Additionally, we find 
that, for each fixed n, the cooperation level in the whole system p c at p = 1.0 is larger 
than that at p — 0. In addition, we would like to point out that it is difficult to use 
theoretical analysis, e.g. pair-approximation method, to investigate the cooperation 
level in this system. However, these simulation results presented in figure 2 clearly 
evidence that the introduction of probabilistic interconnection between two interacting 
networks significantly influences the evolution of cooperation on both of them and there 
exists some intermediate values of p maximizing the cooperation level in the whole 
system. 

In figure 3, we draw the full p — rj phase diagram for the evolution of cooperation on 
the interdependent networks. It is worth noting that, for each value of p, there exists a 
lower critical value and an upper critical value for n respectively. Below the lower critical 
value, defectors dominate the whole population; while above the upper critical value, 
cooperators dominate the whole population. We observe that for intermediate p, both 
the upper and lower boundaries attain their minimum values, where the coupling effect 
between interdependent networks is the strongest. Clearly, the upper critical value of 
n first monotonously decreases until reaching the minimal value at about p = 0.5, then 
increases with increasing p. This result further indicates that the impact of the effective 
interaction topology induced by the probabilistic interconnection reaches the strongest 
at intermediate p. 

In order to intensively investigate the coupling effect induced by probabilistic 
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Figure 4. The frequencies of CC links between two interdependent networks as a 
function of the probabilistic interconnection p with different r\. 



interconnection, we calculate the frequencies of CC links between the interdependent 
networks at equilibrium with different rj in figure 4. Interestingly, we find that the 
frequencies of CC links between the interacting networks with different 77 are pretty 
similar with the frequencies of cooperators in the whole system shown in figure 2(b) 
for corresponding values of 77. This phenomenon indicates that the coupling effect 
due to the internal links between interdependent networks fundamentally influences the 
cooperation level in the whole system. To describe the strategy relationship between 
interdependent networks, especially to describe the assortment of CC links between 
the interacting networks quantitatively, we employ the correlation coefficient between 
individuals [45j 

^E{S'f - E{S % f)^E{S* - E{S 3 f) 

where Si and Sj are the strategies of a random pair of connected individuals i and j. 
If % is a C, Si = 1. Otherwise, Si = 0. E(-) is the expected value of corresponding 
strategy. In an infinite system, we approximately consider a simplified correlation 
coefficient r c = (pec ~ P%)/(Pc ~ Pc)> where pec is the fraction of CC links between 
interdependent networks. We calculate that r c > with different n, when < p c < 1. 
It means that there exists strong positive strategy relationship between interdependent 
networks. Thus, the probabilistic interconnection which fundamentally determines the 
evolution of cooperation in spatial public goods game strengthens the coupling effect 
between the interdependent networks. 

To understand this coupling effect intuitively, in the following we present some 
snapshots of the numerical simulation results of both networks A and B at stationary 
state for fixed n = 0.9 and three different values of p, as shown in figure 5. It can be 
observed that, for intermediate p = 0.5, cooperators can dominate defectors in networks 
A and B, respectively. Whereas for small p = or large p = 1.0, cooperators can 
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Figure 5. Snapshots of the typical distributions of cooperators (white) and defectors 
(black) on network A (a-c) and network B (d-f) obtained by rj = 0.9 and different 
values of p. These snapshots are a 50 x 50 portion of the full 100 x 100 lattices, (a) 
P = ( Pac = 0.7024), (b) p = 0.5 ( Pac = 0.9744), (c) p = 1.0 ( Pac = 0.768), (d) p = 
( P bc = 0.736), (e) p = 0.5 (p bc = 0.9704) and (f) p = 1.0 (p bc = 0.7648). 



only coexist with defectors in the long run. These snapshots demonstrate that, at 
intermediate p, the cooperation level in networks A and B is respectively higher than 
ones for other values of p, thus the fraction of cooperators in the whole population 
is highest at these intermediate p values. Moreover, for different values of p, the 
cooperation level in network A is similar to that in network B. And similar spatial 
patterns simultaneously display in networks A and B especially for intermediate p, 
which indicates the strong coupling effect exhibit between these two interdependent 
networks. 

Finally, to further explore the coupling effect induced by probabilistic 
interconnection between networks A and B on the evolution of cooperation, in figure 6 
we study the time evolution of fraction of cooperators in networks A and B respectively 
by assuming that two different initial conditions for networks A and B. Initially, 
cooperators and defectors are randomly distributed on network A, while only defectors 
are distributed on network B. In figure 6(a), we find that the cooperation level for 
different values of p first decreases sharply and then rapidly increases to the equilibrium 
state. The fraction of cooperators in network B monotonously increases with time for 
different values of p, even for very small p = 0.0001. Interestingly, for each value of p, 
the final fraction of cooperators in network B at equilibrium is similar to the value in 
network A. Remarkably, the fractions of cooperators in networks A and B can easily 
reach an agreement for intermediate interconnection probability p = 0.5. Also, the 
final cooperation levels on both networks for intermediate p are higher than that of 
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Figure 6. Simultaneous time evolution of cooperation level on networks A (p ac ) and 
B (pbc) f° r V = 0.9. Initially, cooperators and defectors are randomly distributed on 
network A with equal probability, but the individuals on network B are all defectors. 



the other values of p. These results imply that the strongest coupling effect between 
the two interdependent networks emerges at intermediate interconnection probability, 
which synchronously promotes the evolution of cooperation in the whole system. 

4. Discussion and conclusion 

Let us further discuss the differences between our model and some relevant previous 
works [33 HQl SHI El HE]. In ref. [3S], Chen et al studied the evolutionary Prisoner's 
Dilemma game on a community network model which exhibits scale-free property. They 
showed that the cooperation level decreases with the increment of the average degree and 
reducing inter-community links can promote cooperation when keeping the total links 
unchanged. In ref. [J7], Wang et al analytically studied the evolution of cooperation 
in multilevel public goods games with community structures. They demonstrated 
that cooperation and punishment are more abundant than defection in the case of 
sufficiently large community size and number with different imitation strength between 
communities. In ref. [38] , Vukov and Szabo studied how the cooperation level is affected 
by the number of hierarchical levels and by the temptation to defect. They showed that 
the highest frequency of cooperation can be observed at the top level if the number 
of hierarchical levels is low, and for larger number of hierarchical levels, the highest 
cooperation level occurs in the middle layers. However, as we described above, these 
works are studied in an isolated single network with community or hierarchical structure 
which does not fully investigate the complex nature of the real world composed of 
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interdependent networks [2H1 EH EES EE] • In our study, the interdependent networks are 
composed of two interacting networks and the internal links build a bridge between the 
two networks. They inspires individuals on different networks to interact with each other 
and then influence each other, thereby promoting the evolution of cooperation in the 
whole system and enhancing the coupling effects between the interdependent networks. 
It is worth pointing out that in our work, although the two interconnected lattice are 
linked directly by the probabilistic interconnection, it does not affect the nature of our 
research. Even if the interaction process and learning process are detached, they can 
still influence the evolutionary dynamics of cooperation in the whole system eventually 
due to the coupling effects between the interdependent networks. 

Recently, Gomez-Gardenes et al studied the evolutionary game dynamics on 
multiplex networks [48]. In their work, each individual plays with all the neighbors on 
different networks and obtain the net payoff of all the payoffs collected in each of network 
layers by using a set of strategies. They showed that the resilience of cooperation for 
extremely large values of the temptation to defect is enhanced by the multiplex structure. 
Furthermore, this resilience is intrinsically related to a non-trivial organization of 
cooperation across the network layers, thus providing a new way out for cooperation 
to survive in structured populations. However, the coupling factors among multiplex 
networks are not considered in their work. Whereas in our study, the introduction of 
probabilistic interconnection controlled by parameter p communicates one network to 
another, introducing the coupling effects. Moreover, the interconnection probability 
p indicates the integration of the two interdependent networks. Unlike the previous 
works, the internal links in our work significantly influence the dynamical behavior in 
the whole system by the coupling effects. There resembles an interesting resonancelike 
phenomenon [191 [50] reflected by the optimal cooperation level at the intermediate 
interconnection probability. Also, introducing appropriate probabilistic interconnection 
between interdependent networks can significantly influence the formation process of 
clusters which determines the evolution of cooperation in the whole system to a great 
extent. Interestingly, the patterns of clusters in both networks are very similar even by 
the approximate locations of cooperators and defectors, which presents an even richer 
resonancelike behavior. Furthermore, comparing the results of p = with p = 1, we find 
that the average level of cooperation when p = 1 is always larger than that of p = 0, thus 
enforcing the positive role of probabilistic interconnection between the two interacting 
networks. Therefore, our study further complements the investigation about coupling 
effect between interdependent networks in the framework of evolutionary graph theory, 
and enriches the knowledge of evolutionary dynamics in the PGG. 

In summary, we have studied the evolution of cooperation on two interdependent 
networks coupled by probabilistic interconnection. We have shown that the introduction 
of probabilistic interconnection provides a new way of understanding the emergence 
and maintenance of cooperation among selfish individuals in sizable groups. We find 
that there exists an optimal intermediate region of p maximizing the cooperation 
level in the whole system. Importantly, we clearly evidence that the introduction 
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of probabilistic interconnection between interdependent networks strengthens the 
coupling effect between them. Therefore, the probabilistic interconnection between the 
interacting networks just like a bridge opening the way for the neighboring networks 
to interact with each other. By means of this simple model, we would like to reveal 
the internal mechanisms about how the evolution of cooperation thrives in the real 
world constructed by interdependent networks. Although this model is simple and does 
not include every kind of circumstances existing, we hope this beneficial attempt can 
highlight the way of exploring the internal mechanisms in promoting the evolution of 
cooperation on interdependent networks which is closer to reality. 
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